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ABSTRACT 

We present a comprehensive analysis of the globular cluster (GC) system of the Local 
Group dwarf irregular galaxy NGC 6822. Our study is based on homogeneous optical 
and near-IR photometry, as well as long-slit spectroscopic observations which are 
used to determine new radial velocities for 6 GCs, two of which had no previous 
spectroscopic information. We construct optical-near IR colour-colour diagrams and 
through comparison to simple stellar population models infer that the GCs have old 
ages consistent with being 9 Gyr or older, while their metallicities are in the range 
between —1.6 < [Fe/H] < —0.4. We conduct a kinematic analysis of the GC population 
and find tentative evidence for weak net rotation of the GC system, in the same sense 
as that exhibited by the underlying spheroid. The most likely amplitude of rotation is 
ss 10 km s _1 , approximately half the magnitude of the observed velocity dispersion. 
Finally, we use the GCs to estimate the dynamical mass of NGC 6822 within ~ 11 kpc 
and we formally find it to be in the range between (3 — 4) x 10 9 Mq. This implies 
an overall mass-to-light ratio in the range of ~ 30 — 40 and indicates that NGC 6822 
is highly dark matter dominated. The mass and the corresponding mass-to-light ratio 
estimates are affected by various additional systematic effects due to limitations of the 
data and the model that are not necessary reflected in the formal uncertainties. 

Key words: Local Group - - galaxies: individual (NGC 6822) — galaxies: kinematics 
and dynamics — globular clusters: general 


1 INTRODUCTION 

According to the ACDM cosmological model, massive galax¬ 
ies build up in part through the amalgamation and merger of 
smaller galaxies. The globular cluster (GC) systems of such 
massive galaxies, which are often quite rich, are also thought 
to be assembled, at least partially, from GCs donated by 
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the smaller accreted galaxies (e.g. ICote et al.lll99Sl . l2000h . 
If this is indeed the case, GCs can act as valuable tracers 
of such events, while offering the potential to understand 
the format ion history of t heir h ost galaxy. In their influen¬ 
tial study, ISearle fc Zinnl (Il978l ) uncovered important clues 
about the formation of the Milky Way through its GC sys¬ 
tem, namely that accretion played a significant part in the 
process. Modern studies of the Galactic GC system arrive 


at the same conclusion (e.g. Zinn 19931: Mackev & Gilmore 

2004 

: Mackev & van den Bergh 20051: Marm-Franch et al. 

to 

0 

0 

: Dotter et al. 2010l. 2011; Keller et al. 20121). Indeed 


the growth of our Galaxy and its GC system is still ongo- 
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ing. Observations of the Sagittarius dwarf galaxy, which is 
currently being accreted onto the Milky Way while donating 
at least 5 GCs to the Galactic halo, provide the strongest 
evidence that this is the case (e.g. iBellazzini et al.l 120031 ; 
lLaw fe MaiewskH 120101 ). The work of our group has shown 
that a similar situation is also observed in M31. The halo of 
this galaxy features a number of prominent stellar streams 
with which many GCs are ass ociated, bearing witne s s to an 
ampl e accretion history (e.g. iMackey et alj 2010bl, |2 013al . 
120141 ; iHuxor et al.ll20lil ; IVelianoski et al. l2013bl . 12014 ). 

The study of dwarf irregular galaxies (dlrr) is a par¬ 
ticularly important topic within the context of under¬ 
standing galaxy formation. Being by far t h e most com ¬ 
mon galaxy type at hig h redshift (e.g. IStiavelli et ai1l2004l ; 
Ivan der Wei et aill2() l 11 1 , dlrr systems could have provided a 
significant contribution toward the growth of massive galax¬ 
ies at early times. In addition, in order to use GC as probes 
of galaxy assembly, it is crucial to understand the properties 
of GC systems in dwarf galaxies, and how they relate to the 
host galaxy properties and to the properties of GCs formed 
in situ around a massive galaxy. In this context, dlrr galax¬ 
ies are especially helpful, since in the local universe they are 
often found either in isolation or in groups without a mas¬ 
sive dominating companion, making them excellent systems 
for studying the pristine properties of GC systems. 

Lying in the southern hemisphere with Galactocentric 
coordinates of l = 25.4°, b = —18.4°, NGC 6822 is one of 
the most i ntriguing dwarf galaxies in the Local Group. Dis¬ 
covered bv lBarnardl tl 18841 4. this barred dlrr galaxy features 
a number of peculiar properties that has made it the focus 
of much attention over the years, even though its low Galac¬ 
tic latitude makes observations quite challenging. It has an 
absolute magnitude My = —15.2 and a h alf-light radius of 
~ 0.5 kpc llHodgelll 977l:lHodge et al.l ll99ll ). Located at a dis¬ 
tance of 472 kpc ( Gorski et al.l 20110 . it is the closest dlrr 
galaxy save for the Magellanic Clouds. It does not appear to 
be associated with either the Milky Way or M31, and it has 
no other neighbouring companions. This galaxy contains a 
substantial spheroidal component, that has been traced via 
its constituent red giant branch star s out to a radius of at 
least ~ 5 kpc dBattinelli et al.l 120061 4 . iDemers et al.l (120061 ) 
used carbon stars in the spheroid as kinematic tracers, find¬ 
ing evidence for rotational support. In addition, NGC 6822 
contains large quantities of HI gas forming a disc, the semi¬ 
major axis of which is positioned perpendicular to the stellar 
spher oid (e.g lde Blok fe Walterll200g . 120061 : IWeldrake et al.l 
I 2 OO 3 ). The HI gas is also observed to be rotating, but at right 
angles to the rotation of the spheroid, showing that there 
are at least two different kinematic components present in 
NGC 6822. This is further supported b y a recent study on 
planet ary nebulae and HII regions by iFlores-Duran et al.l 
(I2014T) . showing that the planetary nebulae have kinematics 
similar to the carbon stars, while the HII regions share the 
kinematics of the HI disc. 

In recent years, dedicated searches have substantially 
increased th e number of GCs known around Local G roup 
galaxies (e.g. lHuxor et al.ll2008l : IVelianoski et al.ll2013al ). in¬ 
cluding NGC 6822. Hubble VII, originally discovered by 
iHubbld dl925l4 . was long considered t he only old GC in 
NGC 6822. However, this changed when lHwang et all (|201lh 
discovered 4 new GCs in wide-held CFHT/MegaCam im¬ 
agery covering an area of 3x3 deg around NGC 6822, . Later, 



Figure 1. Positions of the 8 known GCs overlaid on a Digi¬ 
tal Sky Survey image of NGC 6822 . The colours of the points 
correspond to their radial velocities, as measured in the present 
work (see Section 0 and corrected for the systemic motion of 
NGC 682 2. The ellipse shows t he extent of the stellar spheroid 
traced bv lBattinelli et akl (120061 ). with a semi major axis of 36 ar- 
cmin and an elliptic ity of 0.36. Also sh o wn ar e contours of the HI 
distribution map of lde Blok &; Walter] (l200Cf ) ranging from 3.5 - 
42 xlO 20 cm- 2 . 


IHuxor et alj (|2013h re-examined this data and by comple¬ 
menting it with additional archival Subaru and CFHT im¬ 
agery which expanded the observed footprint, uncovered 3 
additional GCs, bringing the total number of members to 
8. Figure [1] shows the positions of the GCs superimposed 
on an image of NGC 6822, on which are also marked the 
contours of t he HI distribution , rang ing between 3.5 - 42 
xlO 20 cm -2 (Ide Blok fc Walter! 1200(1 1. and the known ex- 
ten t of the red giant branch stars in the galaxy spheroid 
dBattinelli et al.l 120061 ). It is notable that the GC system of 
NGC 6822 is quite spatially extended, with the most remote 
cluster, SCI, having a projected radius of nearly 11 kpc. 
Furthermore, the GCs are arranged in a rather linear con¬ 
figuration, with an orientati on aligned quite closely with the 
majo r axis of the spheroid dHwang et al.ll201ll : IHuxor et al.l 
l2013l) . Five of the 8 GCs in NGC 6822 are fairly extended 
in nature and thus resembl e the extended clusters found 
in the haloes of M31 (e.g . IHuxor et al.l l2005l . 120141) and 
M33 dStonkute et al.l 120081) . The coordinates, position an¬ 
gles, projected radi i (i? Dro i), GC ty pe and Galactic E(B — V) 
extinction values ISchlegel et allfl998l) of the clusters are 
listed in Table |T| 

In this paper we present a comprehensive analysis of the 
GC system in NGC 6822, based on uniform optical and near 
IR photometry and long-slit spectroscopy. The photometry 
is used to constrain the ages and metallicities. Even t hough 
deep e r optical data exist for some clusters (e.g. lHwang et al .1 
l201ll : IHuxor et al] l2013h . the uniformity of the data pre¬ 
sented in this study ensures no systematic offsets are present 
in the age and metallicity estimates due to heterogeneous 
imaging. The spectroscopy yields new radial velocity mea¬ 
surements for 6 of the GCs, two of which had no previous ve- 
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The Globular Cluster System of NGC 6822 3 


Table 1 . Coordinates, position angles, projected radii, GC type (compact or extended), the half light radius (if known) and colour excess 
values due to Galactic interstellar reddening of the GCs in NGC 6822. The adopted centre of NGC 6822 is RA = 19^ 44 m 57.7 s , Dec = 
-14° 48' 12". 


ID 

RA (J2000) 

[h m s] 

Dec (J2000) 

[d m s] 

-^proi 

[kpc] 

9 

[deg] 

GC type 

r h 

[P c ] 

E(B - V) 
[mag] 

Reference 

Hubble-VII 

19 44 55.8 

-14 48 56.2 

0.1 

227 

C 

2.5T0.1 

0.24 

Hubble 11925) 

SCI 

19 40 11.9 

-15 21 46.6 

10.7 

244 

E 

14.0T0.2 

0.16 

Hwanej et al. 12011) 

SC2 

19 43 04.5 

-14 58 21.4 

4.1 

250 

E 

11.5T0.2 

0.22 

Hwan^^et^l. 12011) 

SC3 

19 45 40.2 

-14 49 25.8 

1.4 

94 

E 

7.5T0.5 

0.19 

Hwang et al. 12011) 

SC4 

19 47 30.4 

-14 26 49.3 

6.0 

60 

E 

13.8T0.3 

0.19 

Hwane; et al. 12011) 

SC5 

19 43 42.3 

-14 41 59.7 

2.7 

290 

E 


0.22 

Huxor et al. 12013) 

SC6 

19 45 37.0 

-14 41 10.8 

1.6 

52 

C 


0.19 

Huxor et al. 12013) 

SC7 

19 46 00.7 

-14 32 35.0 

3.0 

43 

C 


0.21 

Huxor et al. 12013) 


locity information. These data are also used to redetermine 
the dynamical mass and mass-to-light ratio of this galaxy. 

2 THE DATA 

2.1 Optical data 

The optical imagery has been obtained as part of the 
PanSTARRSl 3n Survey (Kaiser et al. 2010; K. Cham¬ 
bers et al., in preparation). This survey has targeted three 
quarters of the sky (5 > —30°) in five optical bands, 
gpi, rpi, ipi, zpi, ypi, with the 1.8 m P anSTA RRSl tele¬ 
scope located on Haleakala, Hawaii ( iTonrv et al.ll20lH . The 
optic al system of PanSTARRSl features a 1.4 G igapixel im¬ 
ager (lOnaka et al .1 [20081 : ITonrv fe Oiiakall2QOil) with a field 
of view of 7 deg 2 . The sky has been observed up to four 
times a year in each band, with individual exposures be¬ 
ing between 30 and 45s in length. The data are automat¬ 
ically processe d in real time with the Image Processing 
Pipeline (ipp . lMagnieij|2006l) . After the standard reduction 
processes such as flat-fielding, the pipeline resamples all im¬ 
ages to a uniform pixel size of of 0.25 arcsec and aligns them 
to the equatorial axes onto “skycells” - regular patches on 
the sky of 6250 pixels across. The data are calibrated to each 
other self-consistentl y using partially overlapping exposures 
(ISchlaflv et al.l 120121') , yielding a calibration precision bet¬ 
ter than 10 mmag as measured against SDSS. For the work 
presented here, we use data from the processing version 2 
(PV2) and 2x2 arcmin cutouts of stacked images of the 
GCs provided by the postage stamp server. 

2.2 Near-IR data 

The near-IR imaging of NGC 6822 was taken as part of 
a project to survey red stellar populations of Local Group 
galaxies. The data were o btained in Octob er 2008 using 
the WFCAM instrument (ICasali et al.l 120071) mounted on 
UKIRT. WFCAM has four individual detectors, each with 
a pixel scale of 0.4 arcsec, arranged such that four dithered 
pointings cover a square of 0.75 deg 2 when properly aligned. 
Our observations were conducted in three near-IR bands 
J, H and K. The data were processed with a pipeline cre¬ 
ated by the Cambridge Astronomy Survey UnitQ, which car- 

1 http://casu.ast.cam.ac.uk 


ries out the standard dark current correction, flat-fielding, 
crosstalk removal, sy stematic no i se an d sk y removal. The 
pipeli ne, described in lCioni et al.l (120081 ) and lHodgkin et ahl 
(120091) . also performs full astrometric and photometric cal¬ 
ibration, based on the 2MASS point source catalogue. The 
reduced images are stacked and microstepped to produce 
individual science frames. The finalised frames are not re¬ 
sampled, and thus they retain the original pixel scale of the 
detector. The nightly zero-point variation is < 1%. Figure [5] 
shows 60 x 60 arcsec wide cutouts of each GC hosted by 
NGC 6822 in the PanSTARRSl g P i and UK1RT/WFCAM 
K bands. Note that SCI falls outside the near-IR cover¬ 
age and thus no thumbnail is shown, while SC5 is only 
marginally detected in both optical and near-IR data sets. 

2.3 Spectroscopic data 

Moderate resolution spectra for 6 GCs were acquired in 
the summer of 2009 using the ISIS instrument mounted 
on the 4.2-m William Hershel Telescope ( WHT ), the RC- 
spectrograph installed on the 4-m Kitt-Peak National Obser¬ 
vatory (KPNO) telescope, and the DEIMOS spectrograph 
used by the 10-m W. M. Keck telescope. Several exposures 
of each GC were obtained, with individual integration times 
determined according to the brightness of each target. Our 
observing log is shown in Table [2] 

2.3.1 WHT and KPNO data 

The ISIS instrument comprises two individual detectors, 
attached to separate ‘arms’ of the spectrograph, that in¬ 
dependently sample different wavelength ranges. The blue 
arm covers a wavelength range between ~ 3500 — 5100 A 
with a dispersion of 0.4 A pixel -1 , selected via the usage 
of the EE12 camera along with the R600B grating. The re¬ 
solving power is R ~ 1500. Conversely, the red arm uses the 
RED-PLUS detector together with the R600R grating to 
select the wavelength range between ~ 7400 — 9200 A with 
a dispersion of 0.49 A pixel -1 . The resolving power in this 
case is R ~ 2700. A similar set up was used for the KPNO 
observations. The T2KB detector, in conjunction with the 
KPC007 grating were employed in order to select the wave¬ 
length range between ~ 3500 — 6500 A with a dispersion of 
0.139 A pixel -1 and a resolving power of R ~ 1300. The 
spectra obtained with both telescopes were not binned in 
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Figure 2. gp\ and K band images of the GCs in NGC 6822. Each image is 60 X 60 arcsec wide. North is up and east is left. SCI falls 
outside the near-IR survey and hence no A'-band image is shown, while SC5 is only marginally detected in both passbands. 


Table 2. Log of spectroscopic observations for GCs around NGC 6822 


ID 

Date of obs. 

Number of 

exposures 

Integration time 
per exposure [s] 

Telescope 

Principal 

investigator 

Hubble VII 

14/08/2009 

3 

900 

KPNO 

A. Ferguson 

SCI 

16/08/2009 

3 

1800 

KPNO 

A. Ferguson 

SC3 

17/08/2009 

3 

1800 

KPNO 

A. Ferguson 

SC6 

15/08/2009 

3 

900 

KPNO 

A. Ferguson 

SC7 

15/08/2009 

3 

600 

KPNO 

A. Ferguson 

SC2 

16-17/08/2009 

3 

1800 

WHT 

A. Huxor 

SC7 

16/08/2009 

2 

600 

WHT 

A. Huxor 

SCI 

26/06/2009 

3 

300 

Keck 

J. Hurley 

SC2 

26/06/2009 

3 

300 

Keck 

J. Hurley 

SC6 

26/06/2009 

3 

300 

Keck 

J. Hurley 

SC7 

26/06/2009 

3 

420 

Keck 

J. Hurley 


either the spatial or the wavelength directions. The typical 
S/N of the spectra is ~ 4 — 20 per A. 

The reduction of the spectroscopic data was carried out 
in IRAI0. The basic reduction steps - bias and overscan sub¬ 
traction, flat-fielding, illumination correction - were done 
with appropriate tasks in the CCD package. Having com¬ 
pleted the initial reduction steps on the two-dimensional 
frames, the apall task in the kpnoslit package was used 
to extract one-dimensional spectra from them. We set the 
extraction aperture to be between 2-4 arcsec. With the same 
task, we interactively selected appropriate background sky 
regions, which we fit with a 2nd order Chebyshev polynomial 
and subtracted. The spectra were then traced using a 3rd 


2 IRAF is distributed by the National Optical Astronomy Obser¬ 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 


order cubic spline function, and extracted using the optimal 
variance weighting option in apall. 

We established wavelength calibration of the one¬ 
dimensional frames via Cu-Ne-Ar and He-Ne-AR lamps for 
WHT and KPNO spectra respectively. Comparison ‘arcs’ 
were acquired before and after each cluster exposure. We 
extracted the arc spectra with an identical strategy as the 
target GCs they were used to calibrate. Using the identify 
task, we pinpointed ~ 90 ISIS blue, ~ 25 ISIS red and 
~ 50 RC lines in the arc spectra and fitted a dispersion so¬ 
lution with a 3rd order cubic spline function. The root-mean- 
square residuals of the fits are 0.05±0.01 A, 0 .02±0.01 A and 
0.08 ±0.01 A for the data obtained with the ISIS blue, ISIS 
red and RC instruments respectively. Since two arc spectra 
were observed for each target exposure, the resulting wave¬ 
length solutions were averaged before being applied to the 
appropriate GC via the dispcor task. To assess the reliability 
of the wavelength calibration, we measured the positions of 
sky emission lines in the sky spectra, which are a by-product 
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of the apall extraction. We found that the wavelength cal¬ 
ibration is accurate to 0.08 A for the WHT and 0.1 A for 
the KPNO data. 

Given that all but one cluster were observed as a se¬ 
ries of separate exposures with each instrument, these were 
stac ke d using an inver se variance weighting technique (e.g. 
IVelianoski et ali l2014h . Finally, we continuum subtracted 
the spectra for the purpose of measuring radial velocities. 
Examples of the fully reduced spectra obtained with each 
instrumental set up are shown on Figure [3] The displayed 
spectra are continuum normalized rather than continuum 
subtracted in order to preserve the relative strengths of the 
absorption lines for better visualisation. 

2.3.2 Keck data 

We used the Deep Extragalactic Imaging Multi-Object 
Spectrograph ( DEIMOS ) in longslit mode to observe four 
NGC 6822 GCs (Swinburne Keck program 2009A_W006D, 
PI Hurley). The longslit mode of DEIMOS is achieved by 
inserting a standard single slit mask into the optical path 
of the instrument. For all observations with this spectro¬ 
graph we used a mask featuring a slit 0.8 arcsec wide. The 
DEIMOS camera comprises of 8 CCD chips arranged in a 
4x2 mosaic pattern, with the long axis of this array oriented 
in the spatial direction and the short axis in the dispersion 
direction. This setup essentially produces two spectra simul¬ 
taneously for a given compact target - one on a ‘blue’ CCD 
and one on a ‘red’ CCD. We selected a central wavelength of 
6040 A via the 1200 1 mm -1 grating, resulting in a disper¬ 
sion of 0.33 A pixel -1 . The resolving power of the spectra 
observed with the red CCD is R ~ 4700. 

To reduce observations of a given GC we first identified 
the relevant CCD chips covered by the blue and red spectra. 
We isolated these and reduced them separately; the remain¬ 
ing 6 CCD chips were discarded, as they contained no signal 
coming from our target. Following this, the standard reduc¬ 
tion steps such as bias and overscan subtraction, and flat- 
fielding, were undertaken with appropriate routines in the 
CCD package in iraf. Inspecting the reduced 2-dimensional 
frames, significant geometrical distortions in the dispersion 
of the target objects are clearly evident, as well as distor¬ 
tions in the sky lines. To rectify these we simply wavelength 
calibrated the 2-dimensional frames prior to extracting a 1- 
dimensional spectrum. 

The wavelength solution is based on a Ne-Ar-Kr-Xe 
lamp, taken during the baseline calibrations of the instru¬ 
ment. For the frames containing the redder part of the spec¬ 
trum, we identified 23 lines in the arc spectrum via the iden¬ 
tify and reidentify tasks. The 2-dimensional wavelength so¬ 
lution was then determined with the fitcoords task located 
in the twodspec package. The root-mean-square residuals 
of the fit were 0.02±0.01 A. We next applied the wavelength 
solution to the appropriate frames via the transform routine 
found in the same package. Inspecting the wavelength cal¬ 
ibrated frames, one sees that the distortions in the spatial 
axis have disappeared, and the skylines are perfectly straight 
in a direction perpendicular to the dispersion axis. This is 
a good indicator that the wavelength calibration is reliable. 
For clarity, the wavelength range of our spectra observed 
with the red chip spanned ~ 6000 — 7250 A. 

Unfortunately, we were unable to wavelength calibrate 


the frames containing the blue part of the spectra. This is 
because the arc spectrum observed at these wavelengths con¬ 
tained only three reliable lines, all of which sit extremely 
close to the edge of the frame. This is a well known calibra¬ 
tion issue for data taken in this specific wavelength range 
(albeit one that is nowadays possible to cicumvent). Even 
though this problem forced us to proceed further using only 
the ‘red’ frames, it did not hinder any of the velocity esti¬ 
mation and kinematic analysis which we present in the later 
sections of this paper. 

From the appropriate 2-dimensional frames we ex¬ 
tracted 1-dimensional spectra and stacked these in the same 
manner as for the 4-m data, with the only difference being 
the width of the extraction aperture, which for the Keck 
data we set to 1.1 arcsec. As before, we used the sky spectra 
associated with the extracted 1-dimensional spectra to as¬ 
sess the accuracy of the wavelength calibration. We find that 
the wavelength calibration is typically accurate to 0.06 A. 
An example of a fully reduced, wavelength calibrated Keck 
spectrum is shown in Figure [3] 


3 PHOTOMETRY 

We determined the total magnitudes of each GC in 
NGC 6822 via aperture photometry in all 8 bands, gpi, r P i, 
ip i, zpi, ypi, J, H, and K. The photometric measurements 
were conducted using the phot task in iraf. Accurate centres 
for the classical, i.e. compact, GCs were established with the 
centroid algorithm within phot. Since the astrometric cali¬ 
bration is identical for the entire PanSTARRSl data set we 
performed the centring in the r P i band imaging, since the 
clusters are best seen there. For the diffuse, extended clus¬ 
ters, we determined the centre via visual inspection. Offsets 
from a nearby bright star were used to accurately redeter¬ 
mine the centre of each cluster through the rest of the data. 

We used circular apertures to sum the total flux coming 
from each cluster. The aperture size for a given cluster was 
chosen such that the aperture enclosed the entire extent of 
the target; these sizes are listed in Table [3] The background 
sky was determined between two concentric circular annuli 
surrounding the photometric aperture. We carefully placed 
these in a way so as to exclude unwanted contaminants such 
as extended background galaxies or foreground Milky Way 
stars. The reported magnitude uncertainties, which also in¬ 
clude the variation in the sky, are formally determined by 
the phot task. 

All our photometric measurements were zero-point 
calibrated and corrected for atmospheric extinction, and 
then subsequently corr ected for foreg r ound reddening using 
E(B — V) values from ISchlegel et alj (Il998l l. We calculated 
the extinction coefficients for the PanS TARRSl bands using 
the relations presented in lTonrv et al.l (|2012l 'l: 

Ag P i/E(B—V) = 3.613 - 0.0972(g P i-j P i) + 0.0100(t/ P i-i P i 


(1) 

Ar P i/E(B—V) = 2.585 - 0.0315(p P i-* P1 ) (2) 

Aipi/E(B—V) = 1.908 - 0.0152(g P i-i P i) (3) 

Az P i/E{B-V) = 1.499 - 0.0023(p P i-i P i) (4) 

Ay P i/E(B—V) = 1.251 - 0.0027(s P i-i P i) (5) 
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Figure 3. Examples of fully reduced, continuum-normalized spectra of GCs in NGC 6822. 


ISchlaflv fe Finkbeineil J20H]) recalibrated the E(B — V) 


from Schlegel et al. 1199a ) and recommend that the 


PanSTARRSl reddening coefficients be further multiplied 
by a factor of 0.88, which we included in the magnitude 
derivation. Finally, Table[3]lists the extinction-corrected to¬ 
tal magnitudes of all GCs in NGC 6822. 

It is important to note that SC5 is at best marginally 
detected in the imaging, and thus the measured values repre¬ 
sent lower limits to its magnitudes in all bands (i.e., upper 
limits to its luminosity). The extended cluster SCI is the 
most remote object in this system and is located outside 
the coverage of the near-IR data. 

To allow for comparison with other studies, we convert 
the photometry measurements from the PanSTARRSl sys¬ 
tem to the classical Johnson/C ousins system using equations 
derived in iTonrv et al.l (120121 ') : 


Table 5. ( V — I )p colou r compariso n between the pres ent study 
to those of lHuxor et al. | (120131) and I Hwang et all <1201 ll) . 


ID (V ~ I )0 (V - I )o (V - I )0 

This work Huxo^^L (^20^3) ^20^^) 


Hubble-VII 

0.87 

± 

0.04 







SCI 

0.82 

± 

0.04 




0.85 

± 

0.03 

SC2 

0.90 

± 

0.04 




0.94 

± 

0.03 

SC3 

0.87 

± 

0.05 




1.31 

± 

0.03 

SC4 

1.05 

± 

0.04 




1.12 

± 

0.03 

SC5 

- 0 . 

76 








SC6 

0.87 

± 

0.03 

0.84 

± 

0.03 




SC7 

1.02 

± 

0.03 

1.05 

± 

0.03 





V — rpi = 0.005 + 0.462( fl - r) P1 + 0.013(2 - r) 2 P1 (6a) 
I-i P1 m -0.366 - 0.136((? - r) P1 - 0.018(2 - r) 2 P1 (6b) 

The conversion to the Johnson V and Cousins I bands adds 
an additional 0.012 mag and 0.017 mag to the uncertainty, 
respectively. The results of the conversions are listed in Ta¬ 
ble Q] This table also lists the absolute magnitudes My 0 , 
calculated assuming each cluster has the same distance of 
472 kpc as the centre of NGC 6822. In Table [5] we compare 
our photometric measurements to thos e in previously pu b- 
lished studies. In their discovery paper. iHuxor et al.l (120131 1 
presented optical photometric for SC6 and SC7, measured 


on CFHT/MegaCam archival images. There is an excellent 
agreement between t he photometry presented in this work 
to that in IHuxor et alj (12013l l. with the measurements be¬ 
ing less than one standard deviation apart. Comparing the 
optic al (V — J) q col our from our study to those published in 
IHwang et al.l (I201lh . one finds mixed results. While there is 
an excellent agreement betwee n the colour measure ments of 
our present study and that of lHwang et al.l (120111 1 for SCI 
and SC2, there is a poor agreement regarding the colour 
values of SC4, and even a strong disagreement in the case 
of SC 3. In addition, the colours presented bv iHwang et al.l 
J2011I) are always redder than those in this study. 


© 2015 RAS, MNRAS 000.HHT51 

































































































The Globular Cluster System of NGC 6822 7 


Table 3. Optical and liear-IR photometry of all GCs hosted by NGC 6822. The superscript (1) denotes that all magnitudes for this 
cluster are lower limits (i.e., upper limits to the luminosity). 


ID 

Aperture 

90 

r o 

*0 

z O 

VO 

Jo 

H 0 

Ko 


[arcsec] 

[mag] 

[mag] 

[mag] 

[mag] 

[mag] 

[mag] 

[mag] 

[mag] 


Hubble VII 

3.7 

15.46 

± 

0.02 

14.99 

± 

0.02 

14.72 

± 

0.02 

14.58 

± 

0.03 

14.47 

± 

0.04 

13.44 

± 

0.01 

12.98 

± 

0.02 

12.86 

± 

0.01 

SCI 

7.0 

16.62 

± 

0.02 

16.19 

± 

0.02 

15.97 

± 

0.02 

15.83 

± 

0.03 

15.72 

± 

0.04 










SC2 

8.0 

17.64 

± 

0.02 

16.98 

± 

0.02 

16.80 

± 

0.02 

16.63 

± 

0.03 

16.57 

± 

0.05 

15.53 

± 

0.03 

15.11 

± 

0.03 

14.86 

± 

0.03 

SC3 

3.0 

18.79 

± 

0.02 

18.35 

± 

0.02 

18.08 

± 

0.02 

17.98 

± 

0.03 

17.77 

± 

0.04 

16.93 

± 

0.04 

16.44 

± 

0.04 

16.32 

± 

0.05 

SC4 

8.0 

17.94 

± 

0.02 

17.37 

± 

0.02 

17.00 

± 

0.02 

16.84 

± 

0.03 

16.78 

± 

0.02 

15.74 

± 

0.04 

15.31 

± 

0.04 

15.11 

± 

0.04 

SC5 1 

7.5 

> 

19. 

3 

> 

18. 

6 

> 

18. 

2 

> 

18. 

2 

> 

17. 

8 

> 

18. 

0 

> 

16. 

6 

> 

16. 

4 

SC6 

4.7 

15.76 

± 

0.02 

15.26 

± 

0.02 

15.03 

± 

0.01 

14.91 

± 

0.03 

14.83 

± 

0.02 

13.89 

± 

0.01 

13.42 

± 

0.02 

13.30 

± 

0.01 

SC7 

5.0 

15.24 

± 

0.02 

14.60 

± 

0.02 

14.29 

± 

0.01 

14.12 

± 

0.03 

13.99 

± 

0.02 

12.93 

± 

0.01 

12.39 

± 

0.02 

12.30 

± 

0.01 


Table 4. Photometry of the NGC 6822 GCs converted to the Johnson/Cousins system using Equation l6al and l6bl taken from lTonrv et ahl 
l )2012h . For convenience, the (V — I )o colours are also listed. The superscript is as in Table l3l 


ID 

Vo 

[mag] 

J 0 

[mag] 

(V - I )0 

[mag] 

M V,0 

[mag] 

Hubble VII 

15.05 ± 0.02 

14.18 ± 0.03 

0.87 ± 0.04 

-8.3 

SCI 

16.30 ± 0.02 

15.48 ± 0.03 

0.82 ± 0.04 

-7.1 

SC2 

17.15 ± 0.03 

16.25 ± 0.03 

0.90 ± 0.04 

-6.2 

SC3 

18.44 ± 0.04 

17.57 ± 0.03 

0.87 ± 0.05 

-4.9 

SC4 

17.52 ± 0.03 

16.47 ± 0.02 

1.05 ± 0.04 

-5.9 

SC5 1 

> 18.8 

> 17.6 

~ 0.76 

> -4.6 

SC6 

15.38 ± 0.02 

14.51 ± 0.02 

0.87 ± 0.03 

-8.0 

SC7 

14.77 ± 0.02 

13.75 ± 0.02 

1.02 ± 0.03 

-8.6 


4 AGES AND METALLICITIES 


The best way of determining metallicities and ages of GCs is 
through the analysis of deep colour-magnitude diagrams of 
their resolved stellar populations. Second to that, high qual¬ 
ity spectroscopic data can also be used to derive the metal 
content and age of GCs. However, not all of our spectra have 
sufficiently high S/N to reliably determine these quantities. 
Broadband colour-colour diagrams using a combination of 
optical-optical and optical-near-IR colours are another way 
to estimate ages and metallicities of GCs, and require only 
the existence of homogeneous well-calibrated photometry 
(e.g. Puzia^tjd. 20021: Hempel et alJl2003l : IVelianoski et all 
l2013al : iMunoz et al.l 12014) . The mechanism behind this is 
that both (g — i) o and (g — K)o have similar sensitivity to the 
stars located near the main sequence turn off point, which 
are most sensitive to age, but ( g — K) o is considerably more 
affected by the temperature of the red giant branch which 
is more dependant on metallicity than age. Hence, compar¬ 
ing where GCs on such a plot sit relative to simple stellar 
population model tracks allows a crude estimate of age and 
metal content of each cluster. 

In Figure[4]we show a (gpi — ipi)o vs (gpi-K)o colour- 
colour diagram constructed with the photometry presented 
in Table [3] Integrated magnitudes, in appropriate filters, 
of simple stellar populations having several discrete ages 
are also overlaid on the same plot . Th ey ar e derived from 
the PARSEC1.1 models dBressan et alj|2012h . and span the 
metallicity range of —2.2 < [Fe/H] < +0.6 dex. Formally, 
the diagram shows that the clusters have old ages spanning 
the range between ~ 5 and 12 Gyrs, but they are all consis¬ 
tent with being 9 Gyrs or ol der. This is in agre ement with 
the spectroscopic study by iHwang et al.l (I2014T) . The only 
exception is the cluster SC5, which being only marginally 
detected at best in both the optical and near-IR imaging, 
does not lie near the model tracks and thus its age and metal¬ 


licity can not be reliably constrained. For completeness, it 
is shown as a filled square in Figure [4] 

We estimate the metallicities of the clusters via two 
different methods. In both methods, we ignore SC5 since 
this cluster does not have reliable photometry. In the first 
method, the metallicity of each cluster is determined directly 
from Figure [4] To do this, we projected the position of each 
GC onto a grid of 12 Gyr isochrones in color-color space and 
determined which track lay closest to the GC. The 12 Gyr 
isochrone set is well sampled in metallicity and thus no in¬ 
terpolation is needed for extracting the correct metallicities. 
The upper and lower uncertainties are determined in the 
exactly the same way, considering the bluest and reddest 
possible combination of colours of each clusters respectively, 
as allowed by the colour uncertainties. The metallicities of 
the clusters determined in this way are presented in the first 
column in Tabic [3] 

A different method of estimating GC metallicities is 
through empirical col our-metallicity relat i ons. H ere we use 
the relation derived bv lKissler-Patig et al.l (j2002l ), calibrated 
using 129 Milky Way and M31 GCs that have E(B — V) < 
0.27, to set additional constraints on the metallicities of the 
NGC 6822 GCs. The relation is valid over a metallicity range 
of —2.3 < [Fe/H] < —0.2 dex, and it has the form: 

[Fe/H] = (V — K) o x 1.82 ± 0.11 - 5.52 ± 0.26 (7) 

with an rms of 0.29 dex. Note that such a linear rela- 
tio n between colour and metallicity is a crude assumption 
(ICantiello fe Blakesleell2007l) . 

The metallicities obtained via Equation [7] are also 
shown in Table [6] It is apparent that there is a modest dis¬ 
crepancy between the estimated metallicities via the two dif¬ 
ferent methods. Although all measurements, except that for 
SC2, formally agree with each other, a clear and nearly con¬ 
stant offset of ~ 0.37 dex is present. This mismatch is most 
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(9„1-K)o 

Figure 4. (gpi — ipi)o vs (gpi — K )o colour-co lour diagram of the GCs in NGC 6822, overlaid on top of the simple stellar population 
isochrones derived from the PARSEC 1.1 models (iBressan et a l. 1 I2 oT 3) . Also shown are [Fe/H] values along the 12 Gyr isochrone. SC5 is 
represented by a square since this cluster is a marginal detection in the imaging. The clusters with reliable photometry are all found to 
have old ages. 


likely due to systematics in the stellar evolutionary models. 
Given how rapidly these models progress, at the moment 
we find the metallicities estimated via the empirical relation 
more reliable. 

Table [6] also li sts the [Fe/H] estima tes from the spec¬ 
troscopic study of IHwang et al.l d2014|j . as well as spec- 
troscopic metallicities for Hubble VII from ea rlier studies 
l|Cohen fc Blakesleelll998l : IChandar et alJl200flh . In general, 
the photometric metallicities de rived here are consid erably 
higher than those presented by IHwang et al.l (120146 . with 
only SC3 being in formal agreement. However, no uniform 
offset is present between these two sets of measurements. 
Furthermore, while our photometric [Fe/ H] value for Hubble 
VII is also higher than those der ived in lCohen fc Blakesleel 
dl 99 id) and IChandar et "all d2000l f . the measurements agree 
within the reported uncertainties. 

The reason for the discrepancy bet ween the results pre - 
sented in this contribution and those bv lHwang et al.l d2014l 'l 
is unclear. It is known that NGC 6822 contains young and 
more metal rich populations of stars, and they could be the 
cause of the observed results if they were contaminating the 
photometric apertures. However, repeating the photometry 
with smaller apertures in order to exclude possible contami¬ 
nants did not change the results in Figure^] In addition, the 
ages are recovered to be old which argues against younger 
stellar populations contaminating the derived total magni- 
tu des. As already note d, the optical (V-I) colours presented 
by iHwang et al.1 d201ll 'l are always redder than those in this 
work. It is clear from examining Figure [4] that neither our 
or their optical colours are consistent with their low spec¬ 
troscopic metallicities, at least when interpreted using the 
PARSEC 1.1 models. A colour-magnitude diagram analysis 
of the resolved stellar populations in these GCs would cer¬ 
tainly help to resolve this issue, and might even indicate the 
reason for the mismatch. 


Another effect to be considered is the internal redden¬ 
ing within NGC 6822 itself. As discussed earlier, this galaxy 
contains a significant amount of gas and dust, which, if not 
properly accounted for could cause the clusters to appear 
redder and thus more metal-rich than they actually are. 
Looking at Figure Q] this seems to be unlikely since many 
GCs are located outside the extent of the HI disc, including 
several which are inferred to have high metallicities. 

5 KINEMATICS 

5.1 Radial velocity measurements 

We derived heliocentric radial velocities from the spectra 
obtained with the 4-m class telescopes via a customized y 2 
minimisation routine, using a number of appropriate radial 
velo city standard stars an d bright M31 GCs as templates 
(see I Veljanoski et al.1 12011'1. The derived velocity from each 
science spectrum is the average of all minimisations be¬ 
tween that spectrum and the applied velocity templates. A 
full description of the minimisation process, as well as the 
radial velocity standa rd stars and clusters is presented in 
IVelianoski et al.l (|2014l f. The \ 2 minimisation method is sim- 
ilar to the more frequ ently used cross-correlation technique 
llTonr v & Davis 19796 . with the advantage that it uses the 
uncertainties in the spectra of both the target and template 
objects, which helps to weed out spurious peaks in the \ 2 
functions and reduces the resulting velocity uncertainties. 

During the course of the Keck observations no radial 
velocity standard stars were observed. However, the bright 
cluster SC7 was observed with both the KPNO and the 
WHT telescopes, essentially providing three independent ve¬ 
locity estimates that all agree with each other very well (one 
of these comes from KPNO, while two come from the red 
and blue arms of the ISIS spectrograph on the WHT). Thus 
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Table 6. [Fe/H] values for the GCs with reliable photometry in NGC 6822. Columns refer to: (a) values obtained directly from Figure 
m by determining the projection of each cluster onto the 1 2 Gyr isochrone; (b ) values obtained via Equation [7] an empirical colour - 
me tallicity relation calibrated on t he M31 GC sys tem (Kisslcr-Patie ct al. 2003); spectroscopic metallicity from (c) Hwang ct al.1 d2014h : 
(d) ICohen fe Blakesleel d 1998lj : (e) IChandar et al] d200ol i . 


ID 

[Fe/H] (a) 

[Fe/H] (b) 

[Fe/H] (c) 

[Fe/H] (d) 

[Fe/H] (e) 


[dex] 

[dex] 

[dex] 

[dex] 

[dex] 

Hubble-VII 

_1 14 + 0-04 

A ' i4 — 0.09 

-1.5 ±0.4 

-2.34 ±0.03 

-1.95 ±0.15 

-2.0 ±0.25 

SCI 



-2.00 ±0.04 



SC2 

H R 4 + 0-04 

U.o4_ 0 06 

-1.3 ±0.4 

-2.53 ±0.06 



SC3 

— 1 34+ 006 
0.15 

-1.7 ±0.4 

-1.52 ±0.06 



SC4 

—D 7 c ;+ 0 - 0 6 
u - / o _ 0.05 

-1.1 ±0.4 

-2.53 ±0.08 



SC5 






SC6 

— 1 40+ uu7 
1 * 4U — 0.04 

-1.7 ±0.3 




SC7 

—0 61 +0 03 
u - D 1 - 0.04 

-1.0 ±0.4 





we use our Keck observation of this cluster as a radial ve¬ 
locity template against which the other Keck spectra were 
matched, and adopt as its velocity the error weighted mean 
of the three independent measurements obtained from the 
4-m data sets. The same \ 2 algorithm was used to estimate 
the radial velocity of the Keck targets. 

From all these data we have up to three independent 
velocity estimates for 6 of the 8 GCs in NGC 6822, two of 
which (SC6 and SC7) possess no previous velocity informa¬ 
tion. Table [7] lists the velocity estimates coming from the 
different datasets and their associated uncertainties. Since 
these all agree within the measurement uncertainties, we 
combine them into a single value by performing an error 
weighted average and use these values for the forthcoming 
analysis. For comparison, Table [7] also lists the velocity in¬ 
formation available in the literature, as well as the systemic 
velocity of NGC 6822. Our measurements are generally in 
good agreement with those from the literature, while pos¬ 
sessing considerably smaller uncertainties. 

Of the remaining two clusters, SC5 is sufficiently faint 
that to date there have been no successful v elocit y measure¬ 
ments made for this object. iHwang et al.l (120141 ) present a 
velocity measurement for SC4 (see Table[7]); however the un¬ 
certainties on their estimate are large. Although we did not 
obtain a longslit spectrum of SC4 during any of our observ¬ 
ing runs, individual stars within this object were targeted 
with DEIMOS in multislit mode during our Keck run. The 
purpose of these observations was to investigate the internal 
kinematics of this rather extended cluster, and this analysis 
will be presented elsewhere (Mackey et al., 2015, in prep.). 
Here we adopt the outcome of a computation to determine 
the most-likely radial velocity of this cluster given individual 
measurements for seven member stars, v = —75 ±3 km s . 


5.2 Bayesian kinematic modelling 

Constraining the kinematic properties of the GCs around 
NGC 6822 is one of the main goals of this paper. Unlike in 
massive galaxies, the GC systems of dwarfs are typically 
populated by only a few members, so it is important to 
use a method which will extract the maximum information 
from the sparse data that a re av ailable. Following the exam¬ 
ple from fVelianoski et alj d2014l) . we work in the Bayesian 
framework, which enables us to construct a kinematic model 


and deduce the most likely values for the quantities of in¬ 
terest. 

Given the known rotation of both the HI disk and the 
spheroid (as traced with carbon stars), we particularly want 
to test whether the GC system of NGC 6822 exhibits any 
measurable degree of rotation. We thus construct a kine¬ 
matic model M which features both rotational and veloc¬ 
ity dispersion components. The rotation componen t is con - 
structed following the prescription from lCote et alj (l200ll) : 

u p (0) = u sys + Asin(6 l — 0o) (8) 

where v p are the measured radial velocities, 9 and 9o are 
the position angles of the GCs and the GC rotation axis 
respectively and A is the rotation amplitude. The systemic 
radial motion of the GC system taken to be the heliocentric 
motion of N GC 6822, is labelled as v sys . As elaborated in 
I Cote et~ahi d200ll ). this simple motion along circular orbits 
assumes that the GC system being studied is spherically 
distributed, that the rotation axis lies in the plane of the 
sky and that the projected angular velocity is a function of 
radius only. 

Looking carefully at Figure [T] and Table [7] it is notice¬ 
able that even though the mean velocity of the GCs, which is 
—59 km s -1 , is similar to the systemic motion of NGC 6822, 
there is a preference for the GCs to have radial velocities 
more negative than the systemic velocity of their host. Thus 
we retain v sys as a free parameter instead of fixing in to the 
known systemic velocity of NGC 6822. 

The velocity dispersion is assumed to possess a Gaus¬ 
sian form, and not to change as a function of projected ra¬ 
dius. It is mathematically described as: 

a 2 = ( Av) 2 + a 2 0 (9) 

where oo is the intrinsic velocity dispersion of the NGC 6822 
GC system while An are the uncertainties in the radial ve¬ 
locity measurements. 

Combining these two components, our kinematic model 
A4 takes the following form: 

I _ ( v i ~ )“ 

p i: M(vi,Avi\vp,a) = — -= e 372 (io) 

V2na 2 

where, as before, v p is the systemic rotation described by 
Equation [8] Vi are the observed heliocentric radial veloci¬ 
ties of the GCs as presented in Table [7] and a is a velocity 
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Table 7. Heliocentric radial velocities for the GCs around NGC 6822 as observed from each instrument. Also shown are the combined 
radial velocities which are then used for the kinematic analysis. Note that SC7 was used as a velocity template for the Keck data set. 
For reference, the systemic velocity of NGC 6822, as well as velocity information available in the literature are also shown. 


ID 

KPNO v 

WHT blue v 

WHT red v 

KECK v 

combined v 

Literature v 

Reference for 


[km s -1 ] 

[km s — 1 ] 

[km s -1 ] 

[km s -1 ] 

[km s -1 ] 

[km s -1 ] 

literature data 

NGC 6822 






-57 

McConnachie (£012) 

Hubble VII 

-68 ± 12 




-68 ± 12 

-52; -65 ± 20 

Cohen Sz Blakeslee (19981; Hwang et al. (^014) 

SCI 

-69 ± 16 



-67 ± 4 

-67 ± 4 

-61 ± 20 

Hwanpj et al. (2014) 

SC2 


-72 ± 12 

-80 ± 7 

-74 ± 5 

-76 ± 4 

-106 ± 31 

Hwang et al.^20^^1 

SC3 

SC4 

-83 ± 14 




-83 ± 14 

-71 ± 17 

-115 ± 58; -75 ± 3 

Hwang et al. (2014) 

Hwanpj et al. (2014); Mackev et al. (in prep.) 

SC5 








SC6 

-16 ± 13 



-5 ± 3 

-6 ± 3 



SC7 

-42 ± 17 

-36 ± 11 

-37 ± 2 

template 

-37 ± 2 




dispersion as modelled by Equation [9] Hence, the likelihood 
function for the above kinematic model is: 

Pm(D |0) = C M (-w, Av, 6>|^4, 6» 0 , v sys , cr 0 ) = ( n ) 

i 

where v, Av, and 6 are the observed properties of the GCs, 
while A, do,v B ys and ao are the model free parameters. We 
assume all priors for the free parameters to be flat. 

Even though our model features four free parameters, 
we systematically calculate the likelihood function via a 
“brute-force” method over a regular grid shaped by the pa¬ 
rameter space of the kinematic model. The rotation ampli¬ 
tude A and the velocity dispersion parameters ao range from 
0 to 180 km s' 1 and from 0 to 120 km s'hespectively, both 
with a resolution of 2 km s' 1 . The axis of rotation, do, is 
searched from 0 to 27r rad with a step size of 0.15 rad, while 
the v S ys ranges from -200 to 45 km s' 1 with a step size of 
3 km s' 1 . 

With this set-up, however, it quickly became apparent 
that the algorithm could not properly converge with the po¬ 
sition angle do allowed to vary as a free parameter. This is 
due to the unusual linear configuration of the GCs and the 
small size of the sample. Thus, in order to investigate the 
general kinematic properties of this GC system, we consider 
three constrained models. In the first scenario, we assume 
that if the GC system is rotating, the rotation axis is per¬ 
pendicular to the GC alignment in the plane of the sky. 
This is a reasonable assumption given that rotating config¬ 
urations usually form flattened ellipsoids or disks. We refer 
to this model as to the “disk model”; in this scenario the ro- 
tatio n axis wou ld lie close to that inferred for the spheroid 
bv lDemers et alj ( ;2006l ). In the second case we fix do so that 
the rotation axis lies parallel to the apparent alignment of 
the GCs. This rotating “cigar model” represents the case 
where the cluster system might share the rotation exhibited 
by the HI disk. Finally we consider a scenario in which the 
GCs have kinematic properties separate from those of the 
spheroid and the HI disk, and exhibit no net rotation signal 
(i.e., v p in Equation m is set to zero and the problem is re¬ 
duced to one with only two f ree parameter s). This outcome 
would match that inferred bv lHwang et al.l (l2014fl from their 
kinematic study of SC1-SC4. 

The results of our analysis are shown in Figure [5] The 
general kinematic picture for this GC system is best de¬ 
scribed in the top left panel, where the radial velocities of 
the GCs, corrected for the systemic motion of NGC 6822, 
are plotted against their position angles. In the same panel 


the best-fit rotation curves are overplotted, as described by 
Equation [Sj from our two scenarios which allow for rotation. 
The top right panel shows the marginalized posterior prob¬ 
ability distribution functions for the amplitude in these two 
rotation-enabled models. From these two panels, one can see 
that the disk model detects a low signature of net rotation, 
where the most likely amplitude is 12 ± 10 km s' 1 with 68% 
confidence. The probability of not having any net rotation 
signal is moderately high according to this model. 

For the cigar model, on the other hand, a strong coher¬ 
ent rotation motion is inferred for the GC system. The most 
likely amplitude is 56±31 km s _1 with 68% confidence. This 
is quite a surprising result, since any rapidly rotating system 
ought naturally to form a flattened ellipsoid or disk struc¬ 
ture. One possibility is that we are observing a dynamically 
young system; however, there is at present no evidence that 
this is true for the NGC 6822 GCs - for example one might 
expect to see a stellar stream or other tidal debris associ¬ 
ated with the clusters. A small amount of experimentation 
revealed that this rotation signal is most likely a false detec¬ 
tion driven by SC3 on the one hand, and SC6 and SC7 on the 
other. These objects have strong deviations from the mean 
systemic velocity in opposite directions, and happen to he 
either side of the rotation axis when it is arranged parallel to 
the GC alignment. A small change in the assumed value of 
do in either direction strongly decreases the rotation signal. 
Because there are only a handful of GCs in the entire sys¬ 
tem, these few points are able to substantially influence the 
outcome of the calculation (and indeed we believe it was this 
numerical instability as do moves across the GC alignment 
that led our original four-parameter analysis to fail). 

The lower left and right panels of Figure [5] show the 
marginalized posterior distributions of the velocity disper¬ 
sion do and the systemic velocity of NGC 6822 v sys - Regard¬ 
less of the assumed model, the most likely value of v sys is 
always found to be between -59 and -60 km s _1 with 68% 
confidence intervals of up to ±13 km s' 1 . This is in ex¬ 
cellent agreement with the lit erature v alu e for the systemic 
velocity of NGC 6822 from lMcConnachI3 (l2012l f. The veloc¬ 
ity dispersion is also well defined and does not show a large 
variation depending on the assumed model. Its most likely 
value is 22+ie, 26t]J, 27l 1 g km s _1 for the cigar model, disk 
model, and no-rotation model, respectively. 
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Figure 5. Top left: heliocentric radial velocity, corrected for the systemic motion of NGC 6822, versus position angle of the GCs on the 
sky. The most-likely rotation curves from our two rotation-enabled models are displayed. The disk model shows evidence for a low degree 
of net rotation, while the cigar model detects a more considerable rotation signal (although this is likely to be non-physical - see text). Top 
right: the marginalized posterior distribution functions for rotation amplitude. The amplitude for the disk model peaks at 12d=10 km s 1 , 
where the uncertainty corresponds to the 68% confidence interval. The amplitude for the cigar model peaks at 56 zb 31 km s —1 . Lower 
left: the posterior functions for cro, with peaks between 22-27 km s —1 depending on the model used. In all three cases the distributions are 
quite similar. Lower right: the posterior functions for u S y S - The functions h ave peaks between - 59 and -60 km s —1 , which is in excellent 
agreement with the systemic velocity of NGC 6822 found in the literature (lMcConnachidl2012lI 


6 THE DYNAMICAL MASS OF NGC 6822 

GCs have a long history of being used as probes for estimat¬ 
ing the dynamical masses of the galaxies they surrou nd (e.g. 
iFederici et al.lll990lJl993l ; lEvans &: Wilkinsonll2000l) . We es¬ 
timate the mass of NGC 6822 using its GC s as dynamical 
trace r s following a common approa ch (e.g. I Woodley et al.1 
l20ld ; IVelianoski et al.l l2013bl I2014T ). Since we considered 
three kinematic scenarios, we will use the results of each to 
estimate t he mass of this dwarf gala xy. Solving the Jeans 
equations (iBinnev fe Tremaine! 119871 ') is one common ap¬ 
proach. Since our rotation-enabled models found a varying 
degree of rotation in NGC 6822’s GC system, the total mass 


of this dwarf galaxy can be thought of as the sum from 
two components: a rotating component, and a non-rotating 
pressure-supported component. 

In the scenarios where a net rotation signature is de¬ 
tected, the rotation component, M r , is simply calculated via: 

/? A 2 

M r = (12) 

(jT 

where A is the rotation amplitude of the GC system, i? max 
is the projected radius of the outermost constituent, and G 
is the gravitational constant. 

The pressure supported mass component is deter- 
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mine d using the ‘Tracer Mass Estimator’ (TME) devised by 
lEvans et al. (|2003i ). The TME is mathematically described 
as: 


C 1 v—-\ 

M = (Vi — Vsys) Ri (13) 

i=l 

where R is the projected radius from the centre of NGC 6822 
for a given GC, v is the radial velocity of that GC, corrected 
for global rotation of the system, while v sy s is the systemic 
radial velocity of NGC 6822 itself. The index i loops over 
each GC in the sample of N GCs that have available radial 
velocities. The constant C is dependent on the shape of the 
underlying gravitational potential, a, assumed to be scale- 
free; the radial distribution of the mass tracers, assumed 
to be a power law with index 7 ; and the anisotropy of the 
system. For an anisotropic system, it takes the form: 

4 ( 0 + 7 ) 4-Q-7 l-(nn/r out ) 3 ~ 7 . 

7T 3-7 l-^in/rout) 4 -^ ' 

For 7 = 3, which is typical for a spheroidal stellar halo, 
C becomes: 


c _ 4(q + 3)(1 —q) log(r ou t/n n ) 

n l-(nn/r-out ) 1_<1 


(15) 


We consider cases where the slope of the underlying 
gravitational field, a, is set to 0 , which assumes an isother¬ 
mal halo potential for NGC 6822, and to 0.55, as for a NFW 
profile dNavarro et al.l 1 19961 : IWatkins et al.l l20ld ) . For n n 
and r out we adopt the smallest and largest projected radii, 
respectively, exhibited by the GCs in our sample. As before, 
we include the velocity measured for SC4 by Mackey et al. 
(in prep.); however, we now exclude Hubble-VII from the 
analysis. If a tracer object having a projected radius of zero 
is used in the TME, it will produce a singularity in Equa¬ 
tion ED Hubble-VII has a projected radius of just 0.13 kpc, 
which causes the constant C to increase anomalously. 

Table [ 8 ] lists estimated masses for NGC 6822 within 
~ 11 kpc, i.e. the projected radius of the most remote 
cluster (SCI). The uncertainties of the pressure-supported 
mass component are calculated via the jackknifing tech¬ 
nique, while the uncertainties in the rotation component are 
only due to the error propagation of the detected amplitude 
in the relevant scenarios. Total massess estimated assuming 
the rotating disk and the non-rotating scenario yield similar 
results, which is unsurprising given the low degree of rota¬ 
tion detected in the disk model. For these cases our mea¬ 
surements yield an estimated mass of 3 — 4 x 10 9 Mq, with 
a corresponding mass-to-light ratio of 30 — 40. These esti ¬ 
mates are, however, smaller than those from iHwang et al.l 
<20141) who found the mass of NGC 6822 within 11 kpc to 
be 7.5+Q'f x 10 9 Mq using four extended clusters only (SC1- 
SC4). The discrepancy in the estimates is likely due to the 
large uncertainties in the velocity measurements, in combi¬ 
nation with their smaller cluster sample. Still, our inferred 
mass-to-light ratio makes NGC 6822 a highly dark matter 
dominated system when c ompared to other dwarf galaxies 
of similar luminosity (e.g. lMafedll998l : lMcConnachidl2012l : 
iKirbv et al.ll20l4 ). 

On the other hand, assuming the rotating cigar scenario 


we estimate the mass of NGC 6822 to be unreasonably high, 
about 4 x 10 10 Mq, as indicated by the associated mass- 
to-light ratio. This adds further confidence that the high 
rotation signal found for this model is indeed a spurious 
detection. 

It is worth acknowledging the caveats regarding the 
choice of a and 7 parameters which feature in the TME. It 
is unclear whether their assumed values are indeed the right 
choices when estimating the mass of NGC 6822, thus adding 
additional uncertainty which is not formally included. Us¬ 
ing the same method and identical assumptions, but only 4 
extended GCs, IHwang et al.l (|2014l ) found the total mass of 
NGC 6822 enclosed within 11 kpc to be 7.51 q ® x 10 9 Mq. 
This most likely indicates that when dealing with a small 
samples such as this one, the inclusion or exclusion of a single 
data point can have an important effect on the results. For 
refere nce, using the kinematics of the HI gas. lWeldrake et al.l 
(2003) found a mass of 3.2 x 10 9 Mq within ~ 5 kpc from 
the centre of NGC 6822. 


7 DISCUSSION 


GCs hosted by dwarf galaxies have been often found to 
have similar colou rs as GCs residing in the haloes of mas- 


sive galaxies (e.g 

. Kundu & Whitmore 200ll: lLarsen et al. 

200ll: Peng et al. 

2006), and our results are in keeping with 


that. Comparing the optical colours of the GCs hosted 
by NGC 6822 with those located in the M31 outer halo, 
it is found that they are all mutually consistent, having 
(V — I)o ~ 0.9. This is also true for the GCs residing in 
the dwarf elliptical galaxies NGC 147 and NGC 185, located 
in the outer halo of M31 (IVelianoski et al]|20 13a). In terms 
of their number, extended GCs dominate the population of 
NGC 6822 and there is a similarity between the mean abso¬ 
lute magnitudes of the extended GCs hosted by NGC 6822 
and by M31, which are found to be My 0 = —5.9 ± 0.7 and 
My 0 = —5.6 ±0.7, respectively. This may either be a funda¬ 
mental property of extended GCs in general, or is a result 
of a bias arising from the faint limit of the survey imag¬ 
ing from which these objects were discovered. Given that 
some of them lie right at the detection limit, it is possible 
that the extended GCs have a larger range of absolute mag¬ 
nitudes than observed at present, but the currently avail¬ 
able data allows only for the bright end of the luminosity 
function to be observed, resulting in the low spread and 
similar Mv 0 values. Nonetheless, the range of luminosities, 
structures, ages and metallicities spanned by the GCs in 
NGC6822 is broadly cons istent with that observed in the 
M31 h al o population (e.g. MackjwetjdJ 2006 L 20071 . l2010al . 
l2013al lbl: lAlves-Brito et all 20091 : Sakari et al.l 20151 ). 


A number of studies have focused on GC syst ems hosted 
by d l rr galaxies outside the Local Group (e.g. Seth et all 


|2004 ISharina et al.ll2005l : IGeorgiev et al.ll2006l . 20081 ). They 


have found that GCs have typical optical colours within the 
range of 0.8 < (V — I )0 < 1.1, consistent with our measure¬ 
ments. In their study, IGeorgiev et alJ d2008l ) searched for 
GCs in 19 nearby (2-8 Mpc) dlrr galaxies using archival HST 
data. The galaxies in their sample are members of dwarf 
galaxy associations only, without a dominant massive galaxy 
nearby. In the final sa mple of GC cand i dates that passed all 
their selection criteria, iGeorgiev et al.l l|2008l ) found an ab- 
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Table 8. Dynamical mass estimates for NGC 6822. See text for details. 


Model 

Ngc 

Ol 

M p [M Q ] 

Mp [Mq] 

M total [Mq] 

M/L[M S /L 0 ] 

Disk 

6 

0 

4 ± 2 x 10 9 

3 ± 6 X 10 s 

4 ± 2 x 10 9 

40 

Disk 

6 

0.55 

3 ± 1 X 10 9 

3 ± 6 X 10 s 

4 ± 1 x 10 9 

40 

cigar 

6 

0 

4± 2 x 10 10 

7 ± 8 X 10 9 

4± 2 X 10 10 

400 

cigar 

6 

0.55 

3 ± 1 X 10 10 

7 ± 8 X 10 9 

4 ± 2 X 10 10 

400 

No rotation 

6 

0 

4 ± 1 X 10 9 


4 ± 1 x 10 9 

40 

No rotation 

6 

0.55 

3 ± 1 X 10 9 


3 ± 1 X 10 9 

30 


sence of objects having (V — I )o < 1 and My 0 > —6, i.e. 
faint and blue GCs. The authors claim that this effect is not 
due to a bias in their GC selection criteria, nor a bias due to 
the depth of the imaging, since they are able to identify clus¬ 
ters down to My 0 — —4. In NGC 6822 however, there are 
2 extended GCs which satisfy t hese criteria. If such clusters 
were present in the sample of lGeorgiev et ahl (120081 ). they 
were most likely missed because of their low surface bright¬ 
ness as well as their extended morphology. 

We have considered three scenarios for the kinematic 
configuration of the GC system of NGC 6822. As indi¬ 
cated in the previous Section, the so-called “disk” and “no¬ 
rotation” models are the most plausible of these. In the disk 
scenario our analysis provides tentative evidence for a weak 
rotation signal of amplitude 12±10 km s -1 , although clearly 
the case of zero rotation is not strongly excluded. It is in¬ 
teresting to note that this putative rotation would be in 
the same sense as that detected f or sta rs in the underlying 
stellar spheroid bv lDemers et al.l d2006l ). and also of a com¬ 
parable magnitude (although we note that the Demers et al. 
measurements do not extend to particularly large galacto- 
centric radii). Similarly, the velocity dispersion we infer for 
the GC system mat ches closely that determined from the 
iDemers et al.l 1120061 ) kinematic s ample, once the gl obal ro¬ 
tation has been accounted for ('see lHwang et al1l20 l ll ). These 
similarities suggest that the GC system of NGC 6822 is plau¬ 
sibly associated with its stellar spheroid, and indeed this 
would be consistent with the clear linear alignment of the 
GCs with the major axis of the spheroid. 

Our third kinematic model - the so-called “cigar” model 
- would imply a large systemic rotation amongst the GC 
population, although this is strongly sensitive to the precise 
orientation of the rotation axis. In this case, the sense of 
rotation would be oppos ite to that seen for the HI gas by 
Ide Blok fe Walter! (120001 1: furthermore, the inferred ampli¬ 
tude would be considerably larger than that seen for the HI 
gas at comparable distances along its major axis to those 
of the GCs. In the event this model was correct, it would 
appear that the GCs are not asscociated with the HI gas 
component of NGC 6822. The GCs could constitute a dy¬ 
namically young system, in which case a search for a stellar 
stream along the major axis of the NGC 6822 spheroid may 
prove fruitful. 


8 SUMMARY 

In this paper we present a uniform optical and near-IR pho¬ 
tometric study of the GC system hosted by NGC 6822. The 
photometric measurements are used to estimate the ages 
and metallicities of the GCs via a colour-colour plot and an 
empirical colour-metallicity relation. All GCs are found to 


have ages consistent with 9 Gyrs or older, in agreement with 
past studies. On the other hand, the clusters are found to 
exhibit a range of metallicities (—1.6 < [Fe/H] < —0.4), at 
odds with previous spectroscopic studies which found them 
to be very metal-poor ([Fe/H]< —2.0). Metallicity measure¬ 
ments from resolved colour-magnitude diagrams would alle¬ 
viate this tension, as would higher S/N spectroscopy. 

We also modelled the kinematics of the GC system using 
new velocities for six members, two of which had no previous 
spectroscopic information. With the caveat at our results are 
based on a small sample, we found tentative evidence for a 
weak rotation signal of amplitude 12 ± 10 km s -1 in the 
case where the rotation axis sits perpendicular to the linear 
arrangement of the GCs on the sky. However, zero rotation 
is not strongly excluded by our analysis. If the rotation is 
real, then it would be in the s ame sense as i nf erred for the 
underlying stellar spheroid bv lDemers et al] (l2006h . and of 
a comparable amplitude. Furthermore, the velocity disper¬ 
sion we determine for the GC system, ao = 261^ km s _1 is 
very sim ilar to t hat exhibited by stars in the spheroid (see 
iHwang et al]|2014h . These similarities suggest to us that the 
GCs are plausibly associated with the spheroidal component 
of NGC 6822. 

Finally, using the GCs as dynamical mass tracers, the 
total mass of NGC 6822 is recalculated to be roughly (3 — 
4) x 10 9 Mq, subject to various systematic effects due the 
small data sample and the implicit assumptions built into 
the TME. The corresponding mass-to-light ratio sits in the 
range ~ 30 — 40, implying that NGC 6822 hosts a substantial 
dark matter component. 
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